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The transport in mesoscopic wires with large applied bias voltage has recently attracted great 
interest by measuring the energy distribution of the electrons at a given point of the wire, in Saclay. 
In the diffusive limit with negligible energy relaxation that shows two sharp steps at the Fermi ener- 
gies of the two contacts. Those steps are, however, broadened due to the energy relaxation which is 
relatively weak if only Coulomb electron-electron interaction is assumed. In some of the experiments 
the broadening is more essential reflecting an anomalous energy relaxation rate proportional to E~ 2 
instead of E~ 3 ^ 2 where E is the energy transfer. Later it has been suggested that such relaxation 
rate can be due to electron-electron interaction mediated by Kondo impurity which is, as it is known 
since a long time, singular in E. In the present paper the magnetic impurity mediated interaction 
is systematically studied in the logarithmic approximation valid above the Kondo temperature. In 
the case of large applied bias voltage Kondo resonances are formed at the steps of the distribution 
function and they are narrowed by increasing the bias. An additional Korringa energy broadening 
occurs for the spins by creating electron-hole pairs in the electron gas out of equilibrium. That 
smears the Kondo resonances, and the renormalized coupling can be replaced by a smooth but 
essentially enhanced average coupling and that enhancement can reach the value 8-10. Thus the 
experimental data can be described by formulas without logarithmic Kondo corrections, but with 
enhanced coupling. In certain regions of large bias, that averaged coupling depends weakly on the 
bias. In those cases the distribution function depends only on the ratio of the electron energy and 
the bias, showing scaling behavior. The impurity concentrations estimated from those experiments 
and other dephasing experiments can be very different. An earlier paper has been devoted to provide 
a possible explanation for that by considering the surface spin anisotropy due to strong spin-orbit 
interaction in the host which can hinder the motion of spins with integer spin value. 

PACS numbers: 73.23.-b,72.15.Qm,71.70.Ej 



I. INTRODUCTION 

In the last few years the study of inelastic electron 
scattering mechanisms in mesoscopic metallic systems 
has attracted considerable interest both experimentally 
and theoretically. The motivation is to identify the 
mechanisms responsible for destroying the quantum co- 
herence of the quasiparticles. In one kind of the ex- 
periments the dephasing time is determined from mea- 
surements of magnetoresistancei2iiL2iii, universal conduc- 
tance fluctuation-., and Aharonov-Bohm oscillation in 
mesoscopic rings in magnetic fields^. In the other kind 
of experiments on short wires with large bias voltage the 
non-equilibrium distribution functions of the electron en- 
ergies are measuredi&S which provides an indirect in- 
formation about the energy relaxation of the electrons. 
The dephasing and energy relaxation are due to inelastic 
scattering processes, thus we can obtain information on 
those from these two kinds of experiments. The experi- 
ments were performed at very low temperatures (below 
about IK) where the electron-phonon scattering is weak, 
thus the main inelastic processes are the electron-electron 
scattering and the scatterings by dynamical defects (e.g. 
magnetic impurities, two- level systems). 

In this paper we investigate theoretically the energy re- 



laxation due to magnetic impurities in thin metallic wires 
with large bias voltage. Our motivation is that there are 
many cases where deviations from the expectations of the 
theory of Coulomb electron-electron interaction by Alt- 
shuler, Aronov and coworkers^ were found. The Saclay 
group 7 analyzing the experimental data suggested a phe- 
nomenological anomalous effective energy relaxation rate 
proportional to E~ 2 where E is the energy transfer. That 
differs from the one due to the Coulomb interaction which 
varies like E~ 3 / 2 . That anomalous relaxation rate was 
attributed to the presence of magnetic impurities cither 
implanted or contained by the starting material as con- 
tamination or dynamical defectsii*i2ii&i2ii&ii&. 

The main goal of our earlier paper— has been to sug- 
gest a resolution of the discrepancy concerning the re- 
quired magnetic impurity concentrations to fit the dif- 
ferent experimental data and in some cases those dis- 
crepancies may reach the level of two orders of magni- 
tude. Therefore, even the role of magnetic impurities was 
questioned in spite of the observed sensitivity on external 
magnetic fielcU 7 -. Having a thin wire the magnetic sur- 
face anisotropy^ has been taken into account in Ref. |^ 
which hinders the motion of the spin and in that way 
reduces the Kondo scattering. The previous theoretical 
work9iiii2ii&i2ii&ii&, however, showed the phenomena is 
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rather complex and it is very hard to separate and iden- 
tify the role of the different effects. That is especially 
valid for the non-crossing approximation which is based 
on extensive numerical worki&i^ and many details are 
hidden. The present paper is devoted to those problems, 
and the role of surface anisotropy has already been dis- 
cussed in Ref. 

Considering the Kondo effect the different tempera- 
ture and energy ranges have very different characters. 
Below the characteristic Kondo energy the Fermi liquid 
behavior is recovered. That range has been intensively 
studied even in the non-equilibrium situation in case of 
wires, but the situation in quantum dots shows strong 
similarities 19 ' 20 . The most difficult range is around the 
Kondo energy, while well above that perturbative ap- 
proaches are appropriate, known as logarithmic approx- 
imation. The main idea is to apply perturbation theory 
and to collect the terms with the highest power of the 
characteristic logarithm in each order. That leading log- 
arithmic approach is well established. It will be shown 
that most of the experiments fall in that range. 

In case of non-equilibrium situation the voltage ap- 
plied on the short wire plays crucial role which has to be 
incorporated into that scheme. 

The applied voltage U occurs in two different ways: 

(i) Reduction of the Kondo temperature: In a short 
wire in the diffusive limit with a weak energy relaxation 
the electron comes from one of the two contacts without 
energy loss. The distribution of those electrons is char- 
acterized by different Fermi energies which are split by 
the applied voltage, ell where e is the electronic charge. 
The Kondo effect arises from those regions where there 
are sharp changes in the energy distribution functions, 
thus at the two Fermi energies. At both places separate 
Kondo resonances develop, but the amplitudes of sharp 
changes in the distribution at both energies are reduced 
compared to the equilibrium case, where the change is 
the sum of those two. Therefore, in the wire the equilib- 
rium Kondo temperature for the step with e.g. the lower 
energy is reduced at position x as^i&Si 
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which is valid for Tk < eU where x is measured in the 
units of the length of the wire L. The voltage U measures 
the energy range where the non-equilibrium electron dis- 
tribution takes place. That reduction in Tk can be very 
essential, e.g. for Tk = 0.3 K, eU = 0.3 meV^ 3K, and 
x = 0.5, T K (x = 0.5) = O.ITk = 0.03 K. Thus for a con- 
siderable voltage and for Kondo temperature in a typical 
range (Tk = 0.1 — IK) the reduced Kondo temperature 
is rather small compared to the energy scale eU playing 
role in the electronic transport. That ensures the appli- 
cability of the logarithmic approach as far as Tk <C eU. 
In other cases, e.g. the non-crossing approximation must 
be applie d 13 ' 15 . 

(ii) Enhanced Korringa relaxation rate: the Korringa 
relaxation 22 describes the impurity spin relaxation due 



to electron-hole creations. The creation of those with al- 
most zero energy can be only in the energy range where 
the sharp single-step (T = 0) electron distribution is 
modified. That can be achieved by raising the tempera- 
ture or applying the voltage. In the latter case the extra 
relaxation rate is 
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at T = where 5 is the spin of the impurity and poJ 
is the dimensionless electron-impurity spin coupling. In 
the equilibrium case the Korringa relaxation rate is not 
part of the logarithmic approach, but eU 3> Tk makes it 
non-negligible. 

Namely, in the next to leading logarithmic approxima- 
tion the real part of the pseudospin particle self-energy 
contains a single log which is even lost in the imaginary 
part. That log is not contained by the simple logarithmic 
approach. The effect of the imaginary part (the inverse 
lifetime of the spin) remains, however, very crucial. For 
spins there are different relaxation rates e.g. the T^ 1 
and T 2 _1 in the ESR, NMR experiments. The actual val- 
ues are consequences of a delicate balance between the 
self-energy and vertex corrections, which can modify the 
prefactors like 5(5 + 1) in Eq. (J2J 2 -. Cancellation of 
that type are very well known in the scaling equation of 
the renormalizcd effective coupling, where the value of 
the spin 5 drops ou1i24. Considering the paramagnetic 
impurity electron-spin resonance that has been discussed 
earlier in detail^ 3 - and it is found that in the T 2 _1 relax- 
ation rate the factor (5+1) drops out. Thus the prefactor 
5(5 + 1) in Eq. @ may result in an essential overesti- 
mation for 5 > 1. A systematic study of that problem 
is not known in the present non-equilibrium situation 25 . 
Therefore, the consequences of the present ambiguity is 
only demonstrated by introducing a multiplying factor A 
of order of unity (1/2 < A < 2) for the relaxation rate 
given by Eq. and the role of A will be discussed. It is 
reasonable to assume, that the factor A imitating those 
cancellations is more important for larger spins 5 > 1. 

In the following the logarithmic approach is applied for 
the non-equilibrium case, and the anomalous Korringa 
relaxation is additionally built in. 

Considering the theoretical methods two ways can be 
followed: (i) The Keldysh Green function technique de- 
voted to the non-equilibriumi^ii^iS which may be com- 
bined with the non-crossing approximation 26 , (ii) Direct 
calculation of the scattering rate by the time ordered 
perturbation theory starting with an arbitrary state and 
taking into account the actual distributions in the oc- 
cupations of the intermediate and final states^Si. Re- 
stricting to logarithmic accuracy the second method is 
very easy to apply, in other cases one should turn to the 
first one. 

The paper is organized as follows. In Section ITU the 
transport phenomenon is imposed. In Section IIIII the 
kernel in transport equation due to the magnetic impu- 
rity mediated electron-electron scattering is determined 
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FIG. 1: The wire of length L. The arrow points to the position 
x of the attached tunnel junction which is measured in units 
of L. 



using the leading logarithmic approximation. The elec- 
tron distribution and exchange coupling renormalization 
is carried out in a self-consistent way. In Section IfVI the 
Korringa relaxation is taken into account in the kernel 
and in the smearing of the renormalized coupling, and 
then it is also determined in a self-consistent way with 
the other quantities. In Section [V] the results are ana- 
lyzed from the point of view of the importance of the 
different ingredients of the theory. In the Conclusion the 
results obtained are summarized also from the point of 
view of Ref. 0, where the simplified method is applied 
for the case where the magnetic surface anisotropyi^ is 
included. 



II. NON-EQUILIBRIUM ELECTRON ENERGY 
DISTRIBUTION IN DIFFUSIVE WIRES 

In the experiments the distribution of electron energies 
in a wire with large bias was measured at different posi- 
tion by attaching extra tunnel junction to the side of the 
wire (see Fig. According to that we consider a wire 
with length L and the distribution function at position 
x in the wire in units of L for energy e and at time t is 
denoted by f(e,x,t). The non-equilibrium distribution 
function in the diffusive limit is determined by the time 
dependent Boltzmann equation 

df{e,x,t) 1 d 2 f(e,x,t) 

q7 m h tcoiiXU}) = U (6) 

at td o z x 

where td = jj is the diffusion time and D is the diffu- 
sion constant. In the followings the stationary solution 
is taken, thus the time dependence is dropped and / is 
assumed to be independent of the spin. The collision in- 
tegral lco«.({/}) due to inelastic scatterings in Eq. © 
can be expressed by the scattering rate W(e,E) of elec- 
trons of energy e with energy transfer E as 

hoiiiif}) = J dE{f{e)[\- f{e-E)]W{e,E) 

- [1 — /(e)]/(e — E)W(e — E, —E)}, (4) 

In absence of inelastic scattering processes I C oii.({f}) = 
0, thus the solution of Eq. © is a double-step distribution 
function^ 

f<®(e,x) = (1 - x)n F {e - ^) + xn F (e + ^) (5) 



where U is the voltage between the ends of the wire. 
The sharp steps are smeared due to the inelastic pro- 
cesses even at very low temperature. Taking into ac- 
count inelastic scattering in W and starting with the so- 
lution Eq. |J5J, the Boltzmann equation can be solved 
self-consistently at least numerically. 

In the numerical calculation searching for the station- 
ary solution the Boltzmann equation is solved by itera- 
tion. In the iteration we start with the solution in absence 
of inelastic processes (Eq. and after at each step 
of the iteration the collision integral is calculated from 
the actual distribution function and scattering rate and 
then the new distribution function satisfying the bound- 
ary conditions is determined. The variables x and e are 
discretized usually in 20-60 points. The number of the 
necessary iterations depends on the approximations used 
in the calculation of the collision integral, but usually the 
final solution is reached by 20-120 iteration. 

For electron-electron interaction the scattering rate W 
can be expressed as 

W(e,E) = J de'K(E,e, e')f{e') [l-f( £ ' + E)] (6) 

where the kernel K(E, e, e') is determined by the specific 
interaction and s' is the energy after the collision. 

The predicted dependence of the kernel K on E in 
case of the Coulomb electron-electron scattering^ is ~ 
£;-3/2 m yd regime which explains perfectly some of the 
experiments^, but in other cases other relaxation mech- 
anisms must exist. The Saclay groupi introduced phe- 
nomenologically an extra interaction between the elec- 
trons with K ~ 1/E 2 singularity. It is known since a 
long tirne^S that the interaction mediated by magnetic 
impurities is singular in the energy transfer and recently 
Kaminski and Glazman 1 ^ called the attention to simi- 
lar 1/E 2 singularity in the magnetic impurity mediated 
electron-electron interaction kernel in the leading order 
of perturbation in the Kondo coupling. At that time au- 
thors of Ref. j X ill 1 2ll X Hll X -411 X fill 1 suggested that the energy 
exchange is mediated by Kondo impurities (magnetic and 
structural defects). 

In the following we perform a systematic study of the 
Boltzmann equation for electron-electron interaction 
mediated by Kondo impurities and the different approx- 
imations are built in step by step. 

Assuming low concentration of the impurities we will 
use the single impurity Kondo model— 

H = ^eka^akcr + Jo^Saj^er^/ak/^ (7) 

k.cr kkj 

where a kcr creates a conduction electron with momentum 
k, spin a and energy measured from the Fermi level. 
The conduction electron band is taken with constant en- 
ergy density po for one spin direction, with a sharp and 
symmetric bandwidth cutoff Dq, <j stands for the Pauli 
matrices, Jo is the Kondo coupling, and S is the spin 
operator of the impurity. 
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FIG. 2: The diagrams used for calculating (a) the kernel and 
(b) the Korringa lifetime of the impurity spin. The solid lines 
denote the conduction electrons, the dotted lines the impurity 
spin, and the blob is the Kondo coupling. 



From perturbation theory in the coupling Jo the first 
non-vanishing contribution to the kernel comes from the 
diagram in Fig. [2Ja) asii^ 



K(E,s,s') 
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S(S + l)(2p Jo) 4 
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where c is the concentration of the homogeneously dis- 
tributed magnetic impurities in the wire. That ker- 
nel gives back the phenomenologically introduced energy 
dependence^, but its strength may be too small^ accord- 
ing to the small bare Kondo coupling Jq. 



III. LEADING LOGARITHMIC 
APPROXIMATION 

To improve the calculation as the first step the renor- 
malization of the coupling in the presence of the applied 
voltage is taken into account. The modified kernel reads 

9irr 

K(E,e,e') = -j±plS(S+l) x 

x | J_ (s(S + 1) [J(e) J(s' + E)- J{e')J{e - E)] 2 



+2J(e)J(e')J(e' + E)J{e - E) 



(V 2 (e) + J 2 {e')\J{e' + E)J{e - E) 



E(e -e'-E) 

-[J 2 (e -E) + J 2 (e' + E)]J{e)J(e') I \ (9) 



where J(e) is the renormalized coupling and the approx- 
imations are made in spirit of logarithmic approximation 
discussed in detail in the Appendix. Considering only the 
most divergent terms proportional to l/E 2 in the limit 
E — > and ignore the -E-dependence in the arguments of 
J in Eq. © we get back the result of Eq. (14) of Ref. 12 
thus the kernel is proportional to J 2 (e)J 2 (e')/E 2 . The 
agreement with the result Eq. (10) of Ref. ]TlJ is also only 
after taking the limit E — > in ,/. In general, however, 
this factorization is not true, the kernel contains mixed 
terms as well. 

In our non-equilibrium wire case the renormalized 
Kondo coupling depends also on the position x in the 
wire. In the leading logarithmic approximation carrying 
out similar resummation as in the equilibrium case, the 



leading logarithmic scaling equation is 



dJ(e, x) 

WW) 



2 Po J 2 {f{e-D)- f(e + D)) 



(10) 



where the original bandwidth cutoff Do is reduced to D. 
In Eq. (|10|l |e| D is assumed, thus at the start the elec- 
tron band can be taken always symmetric to e, which as- 
sumption is not applicable in general. On the right hand 
side the actual distribution functions occur which must 
be determined in a self-consistent way. Assuming a sin- 
gle step like distribution that gives back the well known 
result and with a double step the ones in Ref. Q and 
Ref. [3. 

The renormalized Kondo coupling in the leading loga- 
rithmic approximation is the solution of Eq. (|10|) : 



J(e, x) 



Jo 



l-2 Po Jo J ^(f(e-D')-f(e + D')) 




(11) 

It is important to note that due to the finite voltage 
the width of the renormalized, resonant coupling J(e) is 
essentially reduced compared to the equilibrium values 
(see Fig. s |3 and 0}, thus the logarithmic approach can 
be valid even well below the equilibrium Kondo temper- 
ature Tk for large bias U . Due to the finite voltage sep- 
arate Kondo resonances are formediiiiiiiSiiSi at the two 
Fermi energies corresponding to the two contacts thus 
J(s,x) has two peaks (see Fig.s|3]and0}. The validity 
of the logarithmic approach is restricted by the condition 
p J(s,x) < 1 nearby the two Fermi energies which was 
always checked in the calculation by plotting the actual 
Kondo coupling with respect to the energy. 

In case of e U > Tk the Kondo temperature is reduced 
and the singularity is smeared by the Korringa rate, thus 
the condition p^J < 1 is satisfied in most of the cases 
(see Fig.s Eland 0}. 

In the numerical calculation poJ(e,x) > 1 is replaced 
by poJ(e, x) = 1 as poJ(s, x) > 1 is always an artifact of 
the approximation. 



IV. KORRINGA RELAXATION RATE 

In the leading logarithmic order the self-energy cor- 
rections (wave function renormalization and Korringa 
lifetime 2 ^) are neglected as they appear only in higher 
order. At a finite voltage U, however, that inverse life- 
time proportional to U results in anomalous broadening 
which can be very essential. Therefore, the effect of U is 
taken into account in the spin spectral function ps (e) by 
an energy independent Korringa lifetime tk as 



As(0 = - 



2t k 



h 2 



(12) 



The energy dependence in tk leads to changes only in the 
tail of the spectral function which is independent of U for 
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FIG. 3: The renormalized coupling (top) calculated from the 
sharp double-step distribution function (bottom) at x = 0.5; 
T = 0.01K, T K = 0.15K (poJo = 0.039, D = 5 • 10 4 K): 
(a) U — OmV, tk — oo (b) U = 0.15mV, solid line: tk = oo, 
dotted line: ^ = 0.25K (c) U = 0.3mV, solid line: tk = oo, 

dotted line: ^- = 0.25K. 



FIG. 4: The renormalized coupling (top) calculated from the 
sharp double-step distribution function (bottom) at x = 0.5; 
T = 0.01K, Tk = 0.3K (p J = 0.042, Do = 5 • 10 4 K): (a) 
U = OmV, Tk — oo (b) U — 0.15mV, solid line: tk = oo, 
dotted line: ^ = 0.4K (c) U = 0.3mV, solid line: tk = oo, 

dotted line: t^- = 0.4K. 
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\oj\ > eU and therefore it is neglected in the crude leading 
logarithmic approximation applied. The leading term of 
the Korringa relaxation rate is calculated according to 
the self-energy diagram in Fig. H[b) as 

— *=?E f %S(S + l) J tejtfaxHX-ffaxfiffax) 

(13) 

where in the self-consistent calculation the Korringa 
modified couplings Jk are used. In this way the effect of 
U and the position dependence are taken into account. 

To include the effect of Korringa broadening the cal- 
culation of K (E, e, e') given by Eq. © must be repeated 
using the spectral function Eq. I|12|) resulting in the fol- 
lowing form of the kernel: 

2irc 

K(E,e,e',T K ) = —p 3 S(S+l) x 
2E 2 



(E 2 



K 2 \2 



\s{S + 1)[Jk{s)Jk(e' +E)- JK{e')J K {e - E)\ 



+2J K (e)J K (e')J K (s' + E)J K (e - E) 



{E 2 



K 2 \2 



■ (s(S + 1) [J K {s) Jk (e' + E) + J K {e')J K (s - E)] '' 

-2J K (e)J K (e')J K (e' + E)J K (e - E) 
(5(5 + 1) -l)E(e' -e + E) 



{E 2 + £ T )({e> - £ + E) 2 + ^ T ) 

\JW) + Jlie'VKie' + E)J K {e - E) 

-[J 2 K {e -E) + J 2 K (e' + E)]J K {e)J K {e') 

(5(5 + 1)-!)^ 
{E 2 + ^){{e>-E + E) 2 + ^) ' 

[Jk(z) + Jk(£')]Jk(s' + E)J K {e - E) 
+ [J 2 K (e -E) + J 2 K (e' + E)]J K {e)J K {e') 



(14) 



where for the sake of simplicity we suppressed in the nota- 
tion the spatial dependence of the kernel and the Kondo 
coupling. In Eq. and (|14l) Jk{z) is the solution of 
the leading logarithmic scaling equation modified by the 
finite Korringa lifetime 



dJ K (s) 



i; - 2 PQ J 2 K (e) Ids' ' Ps (e-s'){f{e'-D)-f{e'+D)). 



The solution of Eq. (|15|l is the solution of the scaling 
equation Eq. l|llfl with infinite spin lifetime smeared by 
the spin spectral function as 



J^ie)^ / ds'psie-e^J-^e'). 



(15) 



(16) 



Thus at each step of the iteration solving the Boltz- 
mann equation the smeared, renormalized Kondo cou- 
pling Jk(s,x), the distribution function f{e,x), and the 
Korringa lifetime tk (x) were updated. 

The spectral function given by Eq. (|12fl has further 
effects. In the most of the previous calculations the 
Korringa broadening is not taken into account (see e.g. 
Ref. |25|) in the initial and final states. It is easy to 
show that including those broadenings in the collision 
term can be taken into account as an additional broad- 
ening in the singular part of the kernel. As it has already 
been discussed in the Introduction the vertex corrections 
should reduce that broadening. Thus the broadening is 
intensified in that way, but after the reduction by the 
vertex corrections that cannot be very different from the 
overestimated one given by Eq. (|13fl . At the end the am- 
plitude of the relaxation rate is given by Eq. (JTHJ) and 
that modification is taken into account by a phenomeno- 
logical factor A(1/2<A<2) and its influence on the 
results will be discussed later. 



V. RESULTS AND DISCUSSION OF THE 
APPROXIMATIONS 

In the numerical calculation at each step of the itera- 
tion the position dependent collision integral Eq. (@J is 
calculated in terms of the actual distribution function 
and the kernel defined in Eq. © in the following way: 
first the renormalized coupling J is calculated from the 
actual / smeared (Eq. ifTHjl ) or not (Eq. lfTT)l) by the Ko- 
rringa relaxation rate calculated from the preceeding / 
and J. Then the new Korringa relaxation rate and after 
the kernel in terms of tk and J are calculated according 
to the Eq.s {OH) and ifHjl. At the end of the given it- 
eration step the new distribution function is determined 
from the old distribution function and the collision inte- 
gral. The boundary conditions are satisfied at each step 
of the iteration by construction. 

In the followings the results are presented and the im- 
portance of the different ingredients of the theory are 
analyzed. 



A. Renormalized couplings 

The renormalized couplings are calculated using 
Eq. (|llfl or (|16|) without and with the Korringa relax- 
ation for different bias voltages. In Fig.s 01 and 01 the 
couplings calculated from the sharp double-step distribu- 
tion function given by Eq. JSJ) are shown together with 
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the corresponding distribution functions. Without Kor- 
ringa relaxation the resonances formed in the regions of 
the steps are very pronounced. Their widths can be in- 
terpreted as the renormalizcd Kondo temperature. For 
finite biases they are reduced due to the smaller steps 
with sizes x or 1 — x in the distribution function. The 
larger the bias the narrower the resonances are. At the 
overlap of the resonances the coupling can be essentially 
enhanced. 

In Fig. Etb) and (c) (Fig.QJb) and (c)) just for demon- 
stration a fixed value of the Korringa relaxation — 

0.25K (^- = 0.4K) is used. In reality the effect of the 
Korringa broadening is larger for larger bias as t^ 1 is 
linear in U. As it can be seen in Fig.sOIc) and^fc) with 
U = 0.3mV the peaks are almost completely smeared. 
In that case the coupling can be well approximated by 
a constant value poJ, but it is enhanced compared to 
the bare one, e.g. in Fig. [5f a) and (b) J/ Jo ~ 9 and 
J/ Jo ~ 5.5, respectively. In that case the Korringa 
broadening can be estimated according to Eq. (J2J. For 
S = 1/2, x = 0.5, eU = 0.15meV and p J = 0.35 
(eU = 0.3meV and poJ = 0.23) the relaxation rate is 
2^ = 0.22K = 0.19K). For larger bias the depen- 
dence of the approximating constant coupling on the volt- 
age is negligible (see Fig.EJ. (That is the reason why the 
same values of were chosen in Fig. a) and (b) for 
the different biases eU — 0.15meV and eU — 0.3mcV.) In 
this way the Korringa broadening can exceed the width 
of the resonance reduced by the bias and that can be so 
effective that the resonances disappear and the effect of 
the bias is only in the averaged coupling po J. In this way 
the dependence of the final result on is very weak. 

In Fig.[S]we show that the solution of the self-consistent 
calculation using the leading order renormalized coupling 
(solid line) can be well approximated by a constant cou- 
pling (dotted line). The parameters used in the calcu- 
lation are tb = 1.8ns, c = 4ppm, po — 0.21/(site eV), 
S = 1/2, T = 0.05K, T K = 0.15K. 

Based on that observation, the approximation in which 
the self-consistent calculation is performed by using an 
appropriate constant value for the coupling is applied 
in Ref. [3lJ and Paper II where the surface magnetic 
anisotropy for the impurity spin is included. 

It is important to note, however, that the approximat- 
ing constant coupling depends on the voltage bias and 
on the bare coupling (~ Tr-). The former is illustrated 
in Fig. ©for T K = 0.15K. 

The dependence of the distribution function on the 
impurity concentration is shown in Fig. [7| for (a) U — 
0.15mV and (b) U = 0.3mV using the appropriate con- 
stant couplings obtained from Fig. [5] The concentration 
dependent correction is linear only for lower concentra- 
tions. 
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FIG. 5: The coupling constant (top) and the distribution 
function (bottom) obtained from the self-consistent calcu- 
lation using the leading logarithmic renormalized coupling 
(solid line) and an appropriate constant coupling poJ (dot- 
ted line), (a) U = 0.15mV, p J = 0.35, (b) U = 0.3mV, 
poJ = 0.23. The other parameters are td = 1.8ns, c = 4ppm, 
T K = 0.15K, p = 0.21/(site eV), S = 1/2, T = 0.05K. 
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FIG. 6: The dependence of the averaged, approximating cou- 
pling poJ on the voltage bias for Tk ~ 0.15K. The other 
parameters are to = 1.8ns, c = 4ppm, po = 0.21/(site eV), 
T = 0.05K. 
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FIG. 7: The dependence of the distribution functions on the 
impurity concentration is illustrated using the appropriate 
constant couplings obtained (see Fig. [SJ. (a) U = 0.15mV, 
p J = 0.35, (b) U = 0.3mV, p J = 0.23. The other pa- 
rameters are td = 1.8ns, Tk = 0.15K, po = 0.21/(site eV), 
T = 0.05K. 



B. Distribution function 

The calculations are performed according to the iter- 
ation scheme described above. In case of tk = oo the 
kernel given by Eq. ^ is singular in the energy transfer 
E as E~ 2 . There is, however, the cross term with the 
less singular denominator E(s — e' — E). That term is 
less important because of the weaker singularity. Fur- 
thermore that drops out when the energy dependence of 
the couplings is neglected. 




FIG. 8: Comparison of the solutions of the self-consistent 
calculation with and without the cross terms in the kernel 
Eq. 1141 . (a) the distribution functions with (solid line) and 
without the cross terms (dotted line) (b) the difference be- 
tween the two solutions, td = 1.8ns, c = 4ppm, po = 
0.21/(site eV), S = 1/2, x = 0.5, T = 0.05K, T K = 0.3K, 
U = 0.3mV. 

In the presence of finite Korringa relaxation the situa- 
tion is more complex (see the kernel given by Eq. I|14|l). 
The cross term does not disappear even for constant cou- 
pling, but it has limited importance as it is demonstrated 
by calculating the final distribution function without and 
with the cross term. In Fig. the role of the cross terms 
in the kernel given by Eq. I|14[l is illustrated (a) where 
the solution of the self-consistent calculation is plotted 
with (solid line) and without the cross terms in the ker- 
nel (dotted line) . The parameters used in the calculation 
are Trj = 1.8ns, c = 4ppm, po = 0.21/(site eV), S = 1/2, 
T = 0.05K, T K = 0.3K, U = 0.3mV, and x = 0.5. The 
difference is minor which can be seen from Fig. ffib) . De- 
viation occurs only in the vicinity of the steps. Thus in 
most of the cases neglection of the cross terms is a good 
approximation. 



C. The effect of the ambiguity in the strength of 
the Korringa relaxation 

It has been pointed out earlier that the anomalously 
strong Korringa relaxation in the non-equilibrium case 
is not treated in a consequent way as the vertex correc- 
tions and the broadening of the initial and final states 
are not taken into account properly. In order to judge 
their importance a multiplicative factor A is introduced 
as — ► ^27^ (l ater tne notation 7 = XS(S + 1) is 
used). Two sets of curves are calculated for biases 0.1 mV 
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FIG. 9: The distribution functions are illustrated where the 
Korringa relaxation rates are modified by a multiplying fac- 
tor; 7 = XS(S + 1). (a) U = O.lmV (b) U = 0.3mV. The 
other parameters are td = 2.8ns, c = 5ppm, po = 0.21/(site 
eV), T K = 0.3K, T = 0.05K. 



and 0.3mV (see Fig.|5J|. In both cases 7 is changed in the 
wide interval 0.2 < 7 < 2. The differences are relatively 
very small as 7 exceeds the value 7 ~ 0.4. That is a con- 
sequence of the large smearing in the coupling and even 
more in the distribution function. Thus, the improper 
treatment of the vertex corrections is not expected to 
show up as drastic changes in the results. 



D. Distribution function at points out of the 
middle 



The calculation can be performed at any position in the 
wire and the results are demonstrated in Fig. El In these 
cases the reduction in the widths of the Kondo resonances 
are different at the two steps, namely at the smaller step 
it is more reduced. In the distribution function the height 
of the steps is asymmetric according to the actual value 
of x. 



E. Scaling in terms of the applied voltage bias 

Already in the early experiments^ the scaling with re- 
spect to the applied voltage has been observed, thus the 
distribution function at a given position x was described 
by a single variable function f x [-^j). Since that time the 
validity of the scaling has played a central issue both in 
experimental 9 and theoretical studiesAi*i2iili. 





FIG. 10: The coupling constant (top) and the distribution 
function (bottom) obtained from the self-consistent calcula- 
tion at different positions in the wire. The other parameters 
are U = 0.3mV, td = 1.8ns, c = 4ppm, po = 0.21/(site eV), 
T K = 0.15K, T = 0.05K. 



In the present calculation the distribution function at 
position a; is a function of the energy s and voltage bias 
U and it is a functional of the self-consistently deter- 
mined coupling J(e,eU;x), thus it can be written as 
fjj\(s, eU;x). It has been demonstrated in Section IV Al 
that the coupling J{e, e U; x) is a slowly varying function 
of the energy for a fixed value e U in the relevant en- 
ergy window with width e U, which is the dominant one 
concerning the non-equilibrium Boltzmann equation @ 
case. Thus in case of slowly varying J shown in Fig. s [30] 
andJSJthe coupling J[e, x) can be approximated by an av- 
eraged one, J(e U) which is independent of x as it is also 
illustrated in Fig. |3J In that case the distribution func- 
tion is simplified as f{jy(e, e U; x) => f(-^j,J(eU);x). 
Now, if J only slightly depends on the voltage bias 
U (see Fig. |H] for Tk — 0.15K) then the bias depen- 
dence can be dropped for a limited interval of eU and 
J = const, can be used. Then the distribution function is 
/{j}(e, e U; x) =>■ f(-^j,J;x) which exhibits the scaling 
observed experimentally in several cases and illustrated 
in Fig. El 



VI. CONCLUSION 

In this paper the non-equilibrium transport is stud- 
ied in the presence of magnetic impurities representing 
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FIG. 11: The coupling constant (top) and the distribution 
function (bottom) obtained from the self-consistent calcula- 
tion for different biases, (a) x — 0.25, (b) x = 0.5. The 
other parameters are td = 1.8ns, c = 4ppm, Tk = 0.15K, 
po = 0.21/(site eV), T = 0.05K. 



the 1-channel Kondo (1CK) problem. The results can 
be generalized to the 2CK problem which may be real- 
ized by dynamical structural defects^. In the logarith- 
mic approximation scheme applied in the present paper 
there is no difference in the leading logarithmic order as 
the channel index silently follows the continuous electron 
line. The first differences occur in the next approxima- 



tion as closed electron loops appear where summation is 
required with respect to the channel index, thus the Kor- 
ringa relaxation and the vertex corrections are multiplied 
by the number of channels. As it has been discussed pre- 
viously no drastic effects can be expected. 

There is an other drawback of the application of the 
logarithmic approximation in the leading order as that 
overestimates the Kondo temperature by missing the 
prefactor (2po>/) 1 ^ 2 - In general the contribution of the 
next to leading logarithmic approximation can be essen- 
tial, but only nearby the Kondo temperature. The Kor- 
ringa broadening reduces also that effect by reducing the 
averaged coupling p J. On the other hand if the Kondo 
temperature is comparable with the applied voltage, then 
the effective Kondo temperature is less reduced, thus the 
present method cannot be appliedi***^. 

In the present scheme the kernel of the effective mag- 
netic impurity induced electron-electron interaction is de- 
termined. That kernel is governed by the singularity 
E~ 2 in the energy transfer and the less singular cross 
term results only in minor changes. For determination 
of the distribution function a completely self-consistent 
calculation is carried out where in the framework of the 
"poor man's" version of the rcnormalization group the 
couplings are modified by reducing the bandwidth which 
also modifies the Korringa rate and the electron distri- 
bution function. 

For an arbitrary set of the parameters the distribu- 
tion function can be determined using the method de- 
veloped for the weak coupling limit, where Tk is small 
compared to the applied voltage . The relative impor- 
tance of the different ingredients can be judged as the 
method is applied step by step. In the present formula- 
tion some details could be easier followed, than e.g. in 
the non-crossing approximation which is more powerful 
but many details are buried by the huge machinery. 

Considering the comparison with the experiments an 
essential dilemma has been earlier realized. The required 
impurity concentration to fit the experiments in case of 
spin S = 5/2 magnetic impurities is in good accordance 
with values obtained by the other experiments determin- 
ing the dephasing time. In contrary, the cases with pre- 
sumable integer spins show drastic discrepancies. The 
dephasing requires sometimes less than one hundred of 
the impurity concentration and that was attributed to 
the surface anisotropy^ in Ref. 1.3 lL 

The revised version of the form of the surface 
anisotropy 1 ^ makes the effect on the electron distribu- 
tion much weaker, due to its smaller amplitude at larger 
distances measured from the surface. As in Ref. 0^] it is 
shown that the surface anisotropy has weak effect on the 
shape of the distribution, thus the reduced number of im- 
purities with large anisotropy very likely cannot change 
the distribution in an essential way. 

Concidering the dephasing rate, that effect may, how- 
ever, be stronger, as the relevant energy range of the 
electrons is much smaller, as the applied voltage eU is 
replaced by the temperature (eU > kT). Thus, for more 
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FIG. 12: The two electron and one pseudofermion scatter- 
ing amplitude F. The solid and dotted lines represent the 
electrons and the pseudofermion, respectively. One of the 
electron lines connects continuously the upper corners in the 
amplitude F and the other the lower corners. 



impurities, the anisotropy can dominate over the thermal 
energy, and the dephasing rate can be more reduced. 
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APPENDIX 

The derivation of the result Eq. 10 is outlined 
in the non- vanishing leading logarithmic order. The 
Abrikosov's pseudofermion technique 33 is applied where 
the b a operators create spin states with index a and 
then S = ^b^Sapbp. For simplicity the time ordered 

a/3 

diagrams are used with real time. That technique can 
be applied starting with arbitrary states (e.g. non- 
equilibrium) and the transition amplitude between states 
\i) and (/| is given as 



</|T|z) w = £</|ff 
l 



1 



i5 - Hi 



•Hi 



o 



(Al) 



where u) is the energy of the initial state. 

The scattering rate can be calculated in two different 
ways: (i) The scattering amplitudes T are evaluated using 
Eq. (|A.1J| and then those are inserted in the golden rule, 
(ii) The imaginary part of the scattering amplitude is 
combined with the optical theorem (that method is used 
in Ref.[I|. 

In the first method the direction of the spin lines coin- 
cides with the time evolution as only one spin state can 
be present at a given time. In the second method the 
two electron scattering amplitude is determined like in 
the original work of Solyom and Zawadowski^S and the 
spin line is closed. 



(e-e)~' 



tune 



J(co,) 



(e'-rf 



FIG. 13: The two basic diagrams providing singular terms 
in the energy transfer. The solid and dotted lines represent 
the electrons and the pseudofermion, respectively. The square 
box represents the logarithmic electron-pseudofermion vertex 
correction J(u>) and the dotted line the position of the energy 
denominator. 



In the following the first method is used like in the 
papers by Kaminski and Glazmanii and Goppert and 
Grabert. The two basic diagrams providing singular 
terms in the energy transfer E = e — e can be seen in 
Fig. E| where the logarithmic vertex corrections J(u>) 
with two electron and two pseudofermion legs are repre- 
sented by square boxes. In these amplitudes the energy 
is conserved, thus e + e' = e + e' . The energy denomi- 
nators indicated by dotted lines in Fig. ^| separate the 
two vertex corrections T\ and T2 and they are e — e = E 
and e' — s' = e — e = — E where in the latter the energy 
conservation is used. These are singular in the energy 
transfer E. 

As in case of absence of magnetic field and surface 
anisotropy the spins do not carry energy and the in- 
coming and outgoing electrons are on the energy shell, 
thus the energy variables of the vertices are u)\ = e, 
0J2 = e' = e + e' — e = e' + E in Ti and 0J2 = e', 
uji = i = e — E in T2. In terms of the energy dependent 
vertex J(u>) which is real in the logarithmic approxima- 
tion the contributions are proportional to J(e)J(e' + E) 
and J(e')J(e — E), respectively. 

Applying the golden rule the sum of the contributions 
of these two diagrams must be squared. In order to visu- 
alize the result of the detailed calculation we show the two 
forward scattering diagrams in Fig. 1141 in the spirit of the 
optical theorem where the arrows indicate from where the 
imaginary parts are taken, thus the energy is conserved. 
Considering the second "crossing" diagram in the vertex 
on the right hand side the energy variables of the final 
states i <-> e' are interchanged in comparison to the first 
diagram (see Fig. I14|) . That results also in changing the 
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I t 

FIG. 14: The two forward scattering diagrams coming in the 
golden rule. The solid and dotted lines represent the elec- 
trons and the pseudofermion, respectively. The square boxes 
represents the appropriate vertexes (Ti or 1^) and the arrows 
show from where the imaginary parts are taken. 



energy denominators like e — e = E ^ e — e' — — e 1 + e = 
e - e' - E and e' - i' = -E -»• e' - e = -(e - e' - E) 



and also u) 2 = e' + E — > w 2 = § = £ — E in Ti and 
lui = e - E — > uji = e + e' - e = e' + E in T 2 . 

In the "non-crossing" term the energy denominators 
lead to a 1/E 2 singularity. The vertex corrections cor- 
responding to the terms proportional to rf, Y\ an d 
rir 2 (r 2 ri) are J(e) 2 J(e' + E) 2 , J(e') 2 J(e - E) 2 and 
J(e)J(V + E)J{e')J(e - E) occur in Eq. ©. 

In the "crossing" contribution, however, the singularity 
has the form [E{e—e' — E)^ . The corresponding vertex 
corrections are J(e)J(e' + E)J(e)J(s - E), J(e)J(s' + 
E)J{e')J(e' + E), J(e')J(s — E) J(e) J(e — E) , J(e')J{e- 
E)J(e')J(e' + E) in the terms rf, r^, r 2 r 1; and T% 
The vertex corrections with these energy variables occur 
in Eq. ©, as well. 
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